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Muscles shorten faster against light loads than they do against heavy loads. The hyperbolic equation first used by 
A.V. Hill over seven decades ago to illustrate the relationship between shortening velocity and load is still the pre- 
dominant method used to characterize muscle performance, even though it has been regarded as purely empirical 
and lacking precision in predicting velocities at high and low loads. Popularity of the Hill equation has been sus- 
tained perhaps because of historical reasons, but its simplicity is certainly attractive. The descriptive nature of the 
equation does not diminish its role as a useful tool in our quest to understand animal locomotion and optimal 
design of muscle-powered devices like bicycles. In this Review, an analysis is presented to illustrate the connection 
between the historic Hill equation and the kinetics of myosin cross-bridge cycle based on the latest findings on 
myosin motor interaction with actin filaments within the structural confines of a sarcomere. In light of the new 
data and perspective, some previous studies of force-velocity relations of muscle are revisited to further our under- 
standing of muscle mechanics and the underlying biochemical events, specifically how extracellular and intracel- 
lular environment, protein isoform expression, and posttranslational modification of contractile and regulatory 
proteins change the interaction between myosin and actin that in turn alter muscle force, shortening velocity, and 
the relationship between them. 



Muscle cells convert chemical energy into mechanical 
work. Mechanical performance of a muscle is often as- 
sessed by how fast the muscle shortens against a range 
of external loads. Such an assessment is best summa- 
rized by a plot of muscle force versus shortening veloc- 
ity that can be described mathematically by a hyperbolic 
equation of the form: 

(F + «)(V + 4) =c, (l) 

where F and V are force and velocity, and a, b, and c are 
constants. The equation was introduced by A.V. Hill 
(1938), who also suggested that the mechanics of mus- 
cle contraction is closely linked to the muscle's energy 
metabolism, because in his experiments the same hy- 
perbolic force-velocity relationship could be derived 
from heat measurements, and the constant a was found 
to match closely to an empirically derived thermal con- 
stant of shortening heat, a (Hill, 1938). However, a was 
later found not to be a constant but dependent on 
shortening velocity and load (Hill, 1964). It appears, 
therefore, that the force-velocity behavior of a muscle is 
not an unfiltered manifestation of energetic events oc- 
curring inside the muscle, as Hill originally thought. 

Today, the canonical explanation for the characteris- 
tic force-velocity behavior is based on the kinetics of 
cyclic interaction between myosin cross-bridges and 
actin filaments within the contractile units of a muscle, 
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first proposed by A.F. Huxley (1957), followed by im- 
proved models capable of explaining complex behavior 
of the muscle in transient and steady states (Huxley and 
Simmons, 1971; Eisenberg et al., 1980; Pate and Cooke, 
1989; Slawnych etal., 1994; Piazzesiand Lombardi, 1995; 
Edman etal., 1997; Smith etal., 2008; Mansson, 2010; 
Barclay et al., 2010). Enormous amounts of data from 
muscle experiments accumulated over the past decades 
have served to refine the models. In particular, the data 
provided by Piazzesi et al. (2007) has led to a quantum 
increase in our understanding of the molecular basis of 
force-velocity relations in muscle contraction. Interest- 
ingly, the data also provide justification for the use of a 
hyperbolic equation not as a mere empirical description 
but as a meaningful explanation for force-velocity be- 
havior based on cross-bridge kinetics, and may serve to 
revive the Hill equation after decades of under-apprecia- 
tion as a result of incomplete comprehension. 

Force-velocity relationship 

The relationship between muscle force and shortening 
velocity can be visualized by plotting the velocity of a 
shortening muscle as a function of the load (or force) 
pulling on the muscle. A force-velocity curve is usually ob- 
tained from curve-fitting of force-velocity data. The data 
in turn are usually obtained with a protocol involving 
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isotonic quick releases. (A complete description of 
force-velocity relations in muscle should include nega- 
tive loads and positive loads greater than the maximum 
isometric force [F max ] . In this Review, only the force-ve- 
locity relations within the force range of 0 to F max , where 
the Hill hyperbola is most relevant, are considered.) Be- 
fore an isotonic quick release, a muscle is activated at a 
fixed length. The contraction is therefore isometric. 
During an isotonic quick release, the muscle is suddenly 
released from its isometric force to a lower and constant 
force (i.e., isotonic load). In response to the sudden 
change in load, the muscle shortens in a characteristic 
fashion (Civan and Podolsky, 1966) as illustrated in 
Fig. 1. After an initial period of transient changes in 
velocity, muscle shortening settles to a steady velocity. 
The slope of the steady phase (phase 4) of length 
change (Fig. 1 ) is taken as the shortening velocity of the 
muscle corresponding to the isotonic load. Repeated 
velocity measurements from quick releases with differ- 
ent isotonic loads yield force-velocity data points; the 
hyperbolic Hill equation can then be used to fit the data 
and, often by extrapolation, to obtain the maximum ve- 
locity of shortening (V max ) and maximum force (F max ) . 
An example is shown in Fig. 2 for skeletal muscle; simi- 
lar hyperbolic relations are also seen in cardiac and 
smooth muscles. Although the velocity transient (Fig. 1) 
is thought to originate from myosin cross-bridges, it 
should be pointed out that filament compliance could 
alter the transient, especially if it is coupled with a vis- 
cous component. It is therefore essential to identify the 
steady shortening phase (phase 4), preferably with di- 
rect measurement of sarcomere length (Goldman and 
Simmons, 1984; Seowand Ford, 1992), to accurately de- 
termine the velocity of shortening under an isotonic 
load following the transient. 



Length 




Figure 1 . Time course of length response to a force step dur- 
ing an isotonic quick release. The length response can be divided 
into four phases. Phases 1—3 are responses mainly caused by the 
recoil of the elastic muscle elements and cross-bridge transients. 
Typically, within a few milliseconds, the length response settles to 
steady-state shortening (phase 4) characterized by a linear slope. 



Despite its relatively accurate description of the 
force-velocity relations, the Hill equation historically 
has been regarded as purely empirical and devoid of 
any insight into the molecular mechanism of contrac- 
tion (Abbott and Wilkie, 1953). Although this has not 
prevented widespread use of the equation, it may have 
discouraged exploration of physiological significance 
associated with the equation. 

Deviation from Hill's hyperbola at low and high loads 

Careful measurements of force-velocity properties of 
single skeletal muscle fibers have revealed that at low 
loads (less than ^5% F max ), the measured velocities ex- 
ceed those predicted by the Hill hyperbola. At the ex- 
trapolated zero load, Hill's equation underestimates 
the value of V max by ^6-7% (Edman et al., 1976; Julian 
et al., 1986) . For whole muscle preparations with mixed 
fiber types, the underestimation of V max by the Hill 
hyperbola is found to be much greater (Claflin and 
Faulkner, 1989) , but the bulk of the deviation is likely 
caused by heterogeneity in the intrinsic shortening ve- 
locities of the individual fibers within the whole muscle, 
because when the shortening velocity of the whole mus- 
cle exceeds the maximum shortening velocity of the 
slower fibers within the bundle, the shared load in the 
faster fibers increases because the slower fibers can no 
longer keep up with the overall pace of shortening. For 
a muscle preparation with mixed fibers, the force-ve- 
locity relations are better described by a composite of 
different hyperbolic force-velocity curves (Claflin and 
Faulkner, 1989) . 

The Hill hyperbola also failed to accurately describe the 
force-velocity relationship at loads greater than ^80% 
of isometric force (Edman et al., 1976; Lannergren, 
1978). Studies specifically designed to examine the 
deviation of force-velocity data at high loads from Hill's 
hyperbola have revealed double hyperbolic features 
(Edman, 1988; Edman et al., 1997) and reversal of cur- 
vature (Wang et al., 1994; Devrome and Macintosh, 
2007) . Theoretical models of cross-bridge cycling (Edman 
et al., 1997; Negroni and Lascano, 2008; Mansson, 2010) 
were able to simulate the nonhyperbolic features of the 
force-velocity relationship, suggesting that the devia- 
tion from Hill's hyperbola has its origin in the kinetics 
of interaction between myosin cross-bridges and the 
actin filaments. 

Practical applications of the Hill equation 

Despite the discrepancies at low and high loads, velocity 
measurements within the force range of 0.05-0.8 F max 
are remarkably close to the values predicted by the Hill 
equation. Because the maximum power output of a 
muscle occurs at ^0.3 F max , Hill's equation is often used 
in data interpolation to obtain the value of maximum 
power; an example is shown in Fig. 2 (solid circles and 
dashed line) . It is not uncommon to use values of V max 
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Figure 2. Velocity (solid line and open circles) 
and power (dashed line and closed circles) as 
functions of load in an isotonic contraction. 
Modified from Hill (1938) with permission from 
Proceedings of the Royal Society B: Biological Sciences. 
The curve is calculated from the same equation 
(F+ 14.35) (V+ 1.03) = 87.6 used in the original 
fit to data (circles). Power output = F- V. 



and F max extrapolated from the Hill hyperbola as approx- 
imations for the muscle's maximum shortening velocity 
and maximum isometric force. Changes in the curvature 
of the Hill hyperbola, represented mathematically by the 
Hill constant a over the maximum isometric force (i.e., 
a/Fmax), are often used to signify changes in the intrinsic 
properties of the muscle, although the interpretations 
are often vague and superficial because of the lack of un- 
derstanding of the force-velocity behaviors in terms of 
the kinetics of actomyosin interaction. 

Connection between Hill's equation and 
cross-bridge kinetics 

Studies on the molecular basis of cross-bridge kinetics 
and force-velocity relationship (Barany, 1967; Cooke 
and Bialek, 1979; Cooke and Pate, 1985; Nosek et al., 
1987; Cooke et al., 1988; Dantzig et al., 1992; Rayment 
et al., 1993; Potma et al., 1995; Seow and Ford, 1997; 
Piazzesi et al., 2007; Elangovan et al., 2012; Caremani 
et al., 2013) have provided crucial data that allow us to 
make connections between changes in the force-veloc- 
ity properties and the changes in the kinetics of acto- 
myosin interaction. To illustrate this, we need to rewrite 
the Hill equation in a normalized form and compare it to 
an equation derived from actomyosin kinetics. Because 
maximal velocity (V max ) occurs when externally applied 
force (F) is zero, and because at maximum isometric 
force (F max ) the shortening velocity (V) is zero, it follows 
from Eq. 1 that c = (F max + a)b= (V max + b) a, or a/F max = 
tVV max . Therefore, in the normalized form, a single con- 
stant (K) can be used to represent a/F max and b/V msx : 



K = a/¥ m 



(2) 



Now it can be shown that in the normalized form, (F= F/ 
Fmax, V= V/V max ), the Hill equation (Eq. 1) becomes: 



K(l-V) 
K + V 



Force-velocity properties of a muscle are steady-state 
properties. This is because in determining their rela- 
tionship, force and the corresponding velocity are mea- 
sured after the transient response of the muscle has 
settled (Fig. 1). This point becomes important later in 
the analysis of the kinetic model. 

The simplest cross-bridge kinetic model is a two-state 
model, as shown in Fig. 3 A. The entire cross-bridge pop- 
ulation is assumed to reside in either the detached state 
or the attached state, and the corresponding fractional 
cross-bridge populations (D and A) add up to unity; i.e., 
D + A = 1 . Force-velocity relations can be obtained from 
such models because of the intimate relationship between 
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Figure 3. A two-state kinetic model of actomyosin interaction. 
(A) The steady-state distributions of cross-bridge populations 
in detached (D) and attached (A) states are determined by the 
apparent forward (/«>p) and reverse (gxpp) rate constants. (B) A 
three-state model with a detached state (Si), an attached non- 
force-generating state (S2), and an attached force-generating 
state (S 3 ). <?j denotes transition rates from state i to j. 
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force and the fraction of attached cross-bridges, and be- 
tween velocity and the transition rates between the states 
(Huxley, 1957). 

A linear first-order differential equation is used to de- 
scribe the rate of change of the fraction of cross-bridges 
in each state. The rate of change in a particular state is 
taken as the difference between the rates at which the 
cross-bridges move into and out of the state. That is, 



dD l dt = gAPP ■ A ~ IaPP ■ D 



dA/dt = f^p ■ D - gAPP ■ A, 



(4) 



(5) 



where pApp and gXp P are apparent forward and reverse 
transition rates, respectively. In a steady state, dD/dt and 
dA/dt are zero. However, setting Eqs. 4 and 5 to zero 
would result in an ill-posed condition with no unique so- 
lution for the simultaneous equations. This problem is 
avoided by replacing either Eq. 4 or 5 with the following: 



D + A = l. 

From Eq. 4 (with dD/dl = 0), we have: 

a = d (Iapp_ 
\Sapp 

Substituting A in Eq. 6 with Eq. 7: 



D - 



Sapp 



Sapp + )ai 



Similarly, 



Sapp + Iapp 



(6) 



(7) 



(8) 



(9) 



By designating p as force per attached cross-bridge (or 
motor) , the total force (F) produced by the muscle is: 



F = pA = p 



Sapp + J A 



(10) 



To convert Eq. 10 to a hyperbolic function of velocity, 
three assumptions need to be made: (1) the detachment 
rate is linearly proportional to the shortening velocity, 
i.e., gAPP= kV, where k is a constant of proportionality; (2) 
the attachment rate is independent of shortening veloc- 
ity; and (3) force per bridge declines linearly with short- 
ening velocity, i.e., p = 1 — V. With the above assumptions, 
Eq. 10 becomes 



F = (l-V) 



W + f M 



By defining 



K = fAP P /k, 



(11) 



(12) 



and substituting pApp/k with K in Eq. 11, we obtain 
F= (1 — V)K/(V+K), which is exactly the same as Eq. 3 
(the Hill equation). 

The hyperbolic relationship can be demonstrated in 
two-state and multistate models (Baker and Thomas, 
2000; Landesberg and Sideman, 2000; Chin et al., 
2006) . A more general approach in establishing a con- 
nection between the Hill equation and cyclic actomyo- 
sin interaction is illustrated below using a three-state 
model (Fig. 3 B). Assuming that Si is the detached 
state, S> is the attached state before power stroke (non- 
force-generating) , and S3 is the attached state after power 
stroke (force-generating), simultaneous equations can 
be set up to describe the cross-bridge transitions through 
the states: 



dSi/dt = -(c 12 + c 13 ) ■ Sj + cgj ■ S 2 + c 31 • S s 

dS 2 /dt = +c 12 • S 1 - ( c n + c 23 )-S i + c 32 ■ S 3 , 

where cjj is the rate constant from state i to j, and dSi/ dt 
= 0 in a steady state. In matrix form, 

D S=Y, 



where 



D 



1 

-(c 12 +c 13 ) 



(c 21 + c 23 ) c 3 . 



and S = [5/, S 2 , S 3 ] T , and Y = [1, 0, 0] T . Solving for S 3 
(fraction of cross-bridges in the force-generating state) , 



c 12 c 23 + c 13 c 21 + c 13 c 2: 
\D\ 



(13) 



where the determinant \D\ = C21C32 + C21C31 + C23C31 + C12C32 + 
+ + + C13C21 + C13C23. If we assume that the 
transition from S3 to Si (detachment of cross-bridge by 
ATP after the power stroke) is irreversible, then C13 = 0, 
and Eq. 13 becomes 



[c 12 f 31 + C 12 C 32 + %C 3 i + %C 32 + C 23 C 31 ] + Japp 



(14) 



where pAP? = c 12 c 2S . Set gXp P = [c 12 c S i + Ci 2 c 32 + c n c 31 + c 2 ic 32 
+ C23C31], we obtain 



1a, 



Sapp + f&i 



(15) 



Eq. 15 therefore is identical to Eq. 9, and as illustrated 
above for Eq. 9, Eq. 15 can be converted to a hyperbolic 
equation when the same assumptions regarding velocity 



564 Molecular mechanics of muscle 



dependence of gxpp and force per cross-bridge (p) are 
applied. Note that the terms in gxp P all contain one tran- 
sition rate associated with movement of cross-bridges 
away from the force-generating state S3. To obtain a 
hyperbolic force-velocity relationship in a seven-state 
model with multiple force-generating states, only the 
initial attachment rate (from detached state to the first 
attached state) is set to be velocity independent; all 
other transition rates are linear functions of shortening 
velocity (Chin et al., 2006) and as a whole can be repre- 
sented by gxpp as a function of shortening velocity. 

Now let us examine the validity of the three assump- 
tions needed to convert Eq. 10 to the Hill equation in 
light of the findings from Piazzesi et al. (2007) . The first 
assumption is supported by their findings showing that 
the detachment rate increases linearly from ^15/s at 
isometric condition to ^275/ s at a sliding velocity of 
2,000 nm/s, ^2/3 V max (their results are reproduced in 
Fig. 4 A) . If we choose not to ignore the small detach- 
ment rate under isometric condition (i.e., gi = 15/s), 
Eq. 11 becomes: 

F= ^- v {wvf—\ < 16 ) 

{kV + fAPP+giJ 

From Eq. 16 (and Fig. 4 A) it is apparent that at high 
velocities with relatively small compared with k V, the 
equation is essentially the same as the Hill equation. 
But when V approaches zero, overestimation of force 
by the Hill equation becomes obvious. This may ac- 
count for, at least partially, the deviation of force-ve- 
locity data from the Hill hyperbola in the low velocity 
region (Edman, 1988; Wang et al., 1994; Devrome and 
Macintosh, 2007). 

The second assumption partly agrees with the findings 
of Piazzesi et al. (2007). They found that for sliding ve- 
locities greater than ^20% V max (use 3,000 nm/s for the 
value of V max as assumed by Piazzesi et al., 2007), the at- 
tachment rate (Japp) is virtually constant (Fig. 4 B) . At low 
velocities, /app falls below the constant level. For the third 
assumption, the force per attached bridge (p) found by 
Piazzesi et al. (2007) showed linear decline in the high 
velocity region (Fig. 4 C), but again there is discrepancy 
in the low velocity region: the actual value for p under 
isometric conditions falls below the expected value by 
linear extrapolation from the rest of the data (Fig. 4 C) . 

All three assumptions involved in the transformation 
of the kinetic equation (Eq. 1 1) to the Hill equation are 
therefore not entirely correct. However, the discrepan- 
cies are all concentrated in the low velocity region 
where the Hill hyperbola deviates most from the force- 
velocity data. 

Explanations for the deviations from the Hill hyperbola 

As discussed above, if the assumptions regarding the ap- 
parent rates of cross-bridge attachment and detachment, 



and the force per bridge as functions of shortening ve- 
locity are true, then the Hill equation should describe 
perfecdy the relationship between muscle force and ve- 
locity. The fact that the description is not perfect simply 
means that the assumptions are not perfect, but ironi- 
cally, from the deviations some unique behavior of the 
muscle can be understood in terms of molecular mech- 
anisms. The assumption regarding the apparent rate of 
detachment (gxpp) being directly proportional to the 
shortening velocity is true only when the small detach- 
ment rate at isometric condition (g,) is set to zero (see 
Eq. 16 and Fig. 4 A). The implicit assumption embed- 
ded in the Hill equation that g^pp is zero under isometric 
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Figure 4. The apparent rates of detachment (A), attachment 
(B), force per attached bridge (C), and the number of attached 
bridges (D) as functions of shortening velocity. Redrawn from 
Piazzesi et al. (2007) with permission from Cell. 
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conditions (i.e., gi = 0) is therefore one of the sources 
that contribute to the equation's over estimation of 
force at low velocities. A corollary of this analysis is that 
a muscle with higher energetic efficiency in maintain- 
ing isometric force (i.e., small gi) will have less force de- 
viation from the Hill hyperbola in the high force range, 
provided of course all other factors remain the same. 

The apparent rate of attachment (Japp) is constant 
only in the intermediate to high velocity range; at low 
velocities, the rate is slower and not constant (Piazzesi 
et al., 2007) (Fig. 4 B). This means at low velocities the 
fraction of attached bridges is lower than expected 
from extrapolation of the Hill hyperbola (which as- 
sumes pApp remains constant throughout the entire force 
range). This is another source of discrepancy between 
the measured force at the low velocity region and that 
predicted by the Hill hyperbola. To account for the cur- 
vature reversal at low velocities seen in force-velocity 
plots (Wang et al., 1994; Devrome and Macintosh, 2007) , 
Chin et al. (2006) used a /app(v) function (Fig. 5 A) very 
similar to what Piazzesi et al. (2007) found (Fig. 4 B) and 
were able to reproduce the feature (Fig. 5 B) . Chin et al. 
(2006) were also successful in reproducing the double 
hyperbolic feature (Edman, 1988) with a different /app(v) 
function (Fig. 5, C and D) that is also characterized by 
slower attachment rates at low velocities. The more elab- 
orate feature of this function was not observed by 
Piazzesi et al. (2007), but this could be the result of 
insufficient data points in the low velocity region. 

The deviation of p (force per bridge) from the linear 
function in the low velocity region (Fig. 4 C) is yet 
another source of discrepancy between the measured 



force in the low velocity region and that expected from 
the Hill hyperbola (Edman, 1988; Wang et al., 1994; 
Devrome and Macintosh, 2007). Therefore, it appears 
that there are at least three causes that could explain 
why the Hill equation overestimated both the force and 
velocity values in the low velocity (or high force) region 
in the force-velocity relationship. It is likely that a com- 
bination of and interaction among these three causative 
factors result in the deviation of actual measurements 
from the Hill hyperbola. 

It is not surprising that the mid portion of the Hill 
hyperbola fits the force-velocity data quite well, because 
all of the conditions regarding the velocity dependence 
of p, gAPP, and /app inherent to the Hill equation agree 
with experimental evidence (Piazzesi etal., 2007). 

The reason for the underestimation of shortening 
velocity at the other end of the force-velocity curve 
(i.e., at low loads) (Edman et al., 1976; Julian et al., 
1986) by the Hill equation is less clear. No data on p, 
/app, and gxpp in the very high velocity region clearly 
show deviations that could be used to explain the higher 
than expected velocity at low loads. However, that does 
not mean the deviation of the kinetic functions (p(v) , 
Papp(v), and gAPp(v)) from linearity does not exist in the 
high velocity region. 

Force-velocity performance over the range 
of 0.05-0.8 F max 

Most force-velocity studies of muscle performance have 
been performed in a force range that excludes very low 
and high loads. Force-velocity data from these studies 
are generally well described by the Hill hyperbola and 
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Figure 5. Deviation of the apparent at- 
tachment rate (fop) from constancy (dot- 
ted line) at low velocities (A and C) leads 
to the deviation of the simulated force- 
velocity data (circles) from Hill's hyper- 
bola (solid line) at high loads (B and D) . 
Modified from Chin et al. (2006) with 
permission from Biophysical Journal The 
g\pp and force per bridge (p) used in the 
simulation by Chin et al. (2006) were as- 
sumed to be linear functions of velocity. 
See text for more details. 
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from the discussions above; a hyperbolic force-velocity 
relationship is expected and consistent with the results 
from kinetic studies of actomyosin interaction in mus- 
cle cells within this force range. Because the maximum 
power output of muscles occurs within the mid force 
range, variation of this parameter can now be under- 
stood more clearly in terms of the kinetics of actomyo- 
sin interaction. 

Maximum power and the curvature of the force-velocity 
curve. The shape of the force-velocity curve is often de- 
scribed by its curvature, which in a hyperbolic curve can 
be specified by the Hill constant a normalized to F max 
(i.e., a/F max , or Kin Eqs. 2 and 3). As shown in Fig. 6 A, 
an increase in K (or a/F max ) is associated with a decrease 
in curvature, and as shown in Fig. 6 B, a decrease in 
curvature is associated with an increase in the relative 
power. (Relative power is obtained in a force-velocity 
plot in which force and velocity are normalized to 
their respective maximum values.) Because K = fApp/k 
(Eq. 12), it is clear that the curvature is determined by 
the ratio of the apparent attachment rate (Japp) and the 
slope (k) of the gAPp(v) function. It follows that any in- 
tervention that increases /app and decreases the depen- 
dence of gxpp on shortening velocity will make the 
force-velocity curve less curved, and will increase the 
relative power output of the muscle. The increase in 
K value therefore could be caused by an increase in/^p, 
a decrease in k, or a combination of both. Intuitively, an 
increase in the fraction of attached bridges (Eq. 9) 
could lead to an increase in the power output, because 
the power stroke occurs in between attached states. A 
decrease in k also favors accumulation of bridges in at- 
tached states because of the slower rate of detachment 
of cross-bridges as a muscle shortens. A reduced force- 
velocity curvature therefore is related to an increase in 
the fraction of cross-bridges in the attached force-gen- 
erating states at all velocities. Although the explanation 
is derived from a two-state kinetic model, it applies to 
multistate models as well. Similar insight can be derived 
from a three-state model (Fig. 3 B). In a seven-state 
model with five attached and two detached states (Chin 
et al., 2006) , the attachment rate is also found to be cru- 
cial in determining the force-velocity curvature and 
hence relative power. 

The curvature of the force-velocity relation can also 
be derived from Huxley's two-state model (Huxley, 
1957) in terms of kinetic rate constants. It can be shown 
that a/F max is equivalent to the ratio of (fi + gi)/g2 in 
Huxley's model, where fi is related to the attachment 
rate, gi is related to the detachment rate of cross-bridges 
generating positive force, and g2 is the rate constant for 
detachment of compressively strained cross-bridges. Al- 
though there is resemblance to Eq. 12, (fi + gi)/g2 is 
not equivalent to /app/ k. This is because in this particu- 
lar case, Huxley's definition of force-velocity curvature 



is model specific. In a different model, Huxley and 
colleagues (Slawnych et al., 1994) used very different 
rate functions for cross-bridge attachment and detach- 
ment; the definition of the force-velocity curvature for 
this model would therefore be very different from that 
of the Huxley 1957 model. On the other hand, fp^p/k 
can be used to describe the force-velocity curvature as 
long as the data can be satisfactorily described by the 
Hill equation. 

Fraction of attached bridges as a function of velocity. One 
important finding of Piazzesi et al. (2007) is that the 
decrease in force as shortening velocity increases is par- 
tially caused by a decrease in the number of attached 
bridges. Using their /app and gApp as functions of velocity 
in Eq. 9 and ignoring the low velocity region, one can 
see that the fraction of attached bridges (A) decreases 
in a fashion that can be described by a rectangular hy- 
perbola, because /app is constant (except for the low ve- 
locity region) and g^pp increases approximately linearly 
with velocity. Fig. 4 D shows a rectangular hyperbola 
(solid line) fitted to the data from Piazzesi et al. (2007) . 
Because of the relatively good fit and because Eq. 9 is 
part of the link between the Hill equation and the cross- 
bridge kinetics analysis, this implies that an implicit as- 
sumption embedded in the Hill equation is that the 
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Figure 6. Force-velocity (A) and force-power (B) curves with 
different degrees of curvature (a/F max or lvalues). Note that Eq. 3 
can be rewritten in the form: V= K(l — F)/(K+ F). 
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fraction of attached bridges decreases with velocity, as 
Piazzesi et al. (2007) have found experimentally. 

Force-velocity relations in fast and slow muscle fibers. It is 
well known that V max is related to the intrinsic actin- and 
calcium-activated myosin ATPase activity of the muscle 
(Barany, 1967) , which in turn is largely determined by 
the enzymatic properties of the isoforms of the myosin 
family (Moss et al., 1995; Schiaffino and Reggiani, 
1996) . Myosin isoforms in different animal species have 
evolved so that their intrinsic ATPase activities match 
the power requirements imposed by animal body size 
(Seow and Ford, 1991). Within a single species, muscle 
fibers containing predominantly fast or slow myosin iso- 
forms possess, respectively, higher and lower V max values 
(Seow and Ford, 1991; Pellegrino et al., 2003). It is not 
entirely clear how nature modifies ATPase activity, but 
it is likely that the rate of ADP release from its myosin- 
binding site is the target for modification (Barany, 1967; 
Siemankowski etal., 1985; Raymentetal., 1993; Elangovan 
et al., 2012). This implies that at maximal shortening 
velocity, the rate-limiting step in a cross-bridge cycle is 
the hydrolysis of ATP and the release of P; and ADP 
from the myosin head. The fast and slow myosin iso- 
forms, besides determining V max , also shape the force- 
velocity curve. Woledge et al. (1985), after compiling 
force-velocity data from many independent studies 
(Katz, 1939; Close, 1964; Woledge, 1968; Cecchi et al., 
1978; Lannergren, 1978; Luff, 1981; Lannergren et al., 
1982; Ranatunga, 1982), found a positive correlation 
between a/F max and V max . Later studies examining short- 
ening velocities from fast and slow muscle fibers also 
showed a similar correlation (Brooks and Faulkner, 
1988; Stienen et al., 1988; Barclay et al., 1993; Wahr and 
Metzger, 1998; Gilliver et al., 2009). This indicates that 
a fast muscle like extensor digitorum longus (EDL), 
compared with a slow one like soleus, not only shortens 
faster but also has greater power output because of less 
curvature in its force-velocity curve. According to Eq. 12, 
in order for a fast muscle to have less curvature in its 
force-velocity relation, either the apparent attachment 
rate (Japp) of the myosin to actin needs to be increased 
or the rate of increase in £app as a function of increas- 
ing velocity needs to be decreased, or a combination 
of both. 

In cardiac myocytes containing a-MHC and (3-MHC 
isoforms, it has been found that the faster isoform 
(a-MHC) is associated with higher shortening velocity 
and power output, as well as a greater value for a/F max , 
when compared with the slower (3-MHC isoform 
(Herron etal., 2001). 

Variation in force-velocity properties exists within a 
single fiber type of the same species and also across 
species. Examination of contractile properties of type I 
fibers from different species has revealed the same posi- 
tive correlation between a/F max and V max (Seow and 



Ford, 1991; Widrick etal., 1997). Interestingly, the re- 
lationship between a/F max and V max obtained from a 
single fiber type of the same species is opposite to that 
from different fiber types or the same type that be- 
longs to different species (Gilliver et al., 2009, 2011); 
that is, the ratio a/F max decreases with increasing V max . 
This indicates that modification of actomyosin kinetics 
resulting in variation of force-velocity properties within 
a fiber type of a single species is different from the 
modification that occurs in the fast and slow types in 
the same species, or within the same type from differ- 
ent species. 

Effects of temperature. In summarizing the data from 
Lannergren et al. (1982) and Ranatunga (1982) , Woledge 
et al. (1985) noticed that V max is more sensitive to tem- 
perature than F max is. Because shortening velocity is 
mosdy determined by the apparent detachment rate 
(gXpp) , which in turn is closely related to the rate of ADP 
release after the power stroke under low loads, the high 
temperature sensitivity of V max is likely a reflection of 
the same temperature sensitivity of actomyosin inter- 
action, especially in the conformational change step 
within the myosin molecule leading to ADP release 
(Rayment et al., 1993). Ranatunga (1984) specifically 
examined the temperature effect on V max and a/F max in 
fast (EDL) and slow (soleus) muscles from rats. At each 
of the six temperatures examined, a/F max was found to 
correlate positively with V max , and across a temperature 
range of 10-35°C, the value of a/F max increases with in- 
creasing temperature in both the fast and slow muscles 
(Ranatunga, 1984). The increase in a/F max with tem- 
perature could be a secondary effect of increasing V max 
with temperature, because as discussed earlier, a/F max 
and V max are positively correlated. However, a closer ex- 
amination of the data from Ranatunga (1984) suggests 
that this is not the case. Fig. 7 shows that although the 
values of a/F max in both slow and fast muscles increase 
with temperature, they do not have the same velocity 
dependence. The «/F max value for a fast muscle in- 
creases faster with velocity compared with the relation- 
ship for a slow muscle (Fig. 7) . The increase in a/F max 
with temperature is therefore muscle-type specific. The 
«/F max value of the fast muscle has an overall greater 
sensitivity to temperature compared with that of slow 
muscle. Increasing temperature does not merely in- 
crease the shortening velocity. If the velocity values in a 
force-velocity curve are scaled up by a constant factor, 
the a/F max value does not change. The fact that a/F max 
is altered by temperature suggests that /app, k, or both 
are temperature sensitive, at least in mammals. Zhao 
and Kawai (1994) have found that the fraction of de- 
tached cross-bridges declines as temperature increases, 
indicating an increase in /app (and thus a/F max ) with in- 
creasing temperature and providing an answer to the 
question of why force-velocity curvature decreases with 
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increasing temperature. In some cold-water fish species it 
has been found that a/F max is not temperature sensitive, 
even though V max increases with temperature (Johnston 
and Salamonski, 1984; Johnston and Sidell, 1984). Fish 
myosin isoforms may have adapted to temperature 
ranges to minimize loss of power at low temperatures. 

Effects of[ADP], [ATP], and fPJ. An increase in ADP con- 
centration surrounding the site of actomyosin interac- 
tion will lead to a decrease in gxp P because it interferes 
with the release of ADP from the cross-bridges (Barany, 
1967; Cooke and Pate, 1985; Siemankowski et al., 1985; 
Seow and Ford, 1997; Elangovan et al., 2012). The ex- 
pected effects are an increase in the fraction of attached 
bridges (therefore isometric force) and a decrease in 
shortening velocity. These effects have been observed 
in experiments (Cooke and Pate, 1985; Seow and Ford, 
1997) and simulated with a multistate model of cross- 
bridge kinetics (Chin et al., 2006). There is an increase 
in the a/F max value as a result of the increase in [ADP] 
(Seow and Ford, 1997) . Because reducing the rate of 
ADP release directly affects gxp P , and because of the 
close relationship between gxp P and shortening velocity 
(Piazzesi et al., 2007), the increase in a/F max (or reduc- 
tion of curvature of the force-velocity curve) is likely 
caused by a reduced dependency of gXp P on shortening 
velocity (i.e., a smaller lvalue) . It should be pointed out 
that the ATP-regenerating system (phosphocreatine 
and creatine kinase) is omitted in experiments examin- 
ing the effects of ADP to control ADP concentration. 

A decrease in [ATP] has very similar effects as an in- 
crease in [ADP] without the ATP-regenerating system. 
It causes an increase in isometric force and a decrease 
in shortening velocity (Ferenczi et al., 1984; Cooke and 
Pate, 1985; Stienen et al., 1988; Seow and Ford, 1997), 
and an increase in a/F max (Ferenczi et al., 1984; Stienen 
et al., 1988; Seow and Ford, 1997; Wakayama and 
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Figure 7. The relationship between «/F max and V max of fast 
(EDL) muscle (closed circles) and slow (soleus) muscle (open 
circles) at different temperatures. Data source from Ranatunga 
(1984) with permission from The journal of Physiology. 



Yamada, 2000). An increase in [ADP] associated with 
the experimental conditions creating low [ATP] could 
be an explanation for the observed effect of ATP 
(Cooke and Bialek, 1979) . But even without an increase 
in [ADP] as a secondary effect, lowering [ATP] will re- 
duce gxpp and hence give rise to an effect similar to that 
caused by increased [ADP]. 

Amrinone (a bipyridine compound) has been shown 
to increase ADP affinity to actomyosin (Albet-Torres 
et al., 2009). Not surprisingly, the effects of amrinone 
on the force-velocity relationship are very similar to the 
effects of high [ADP] described above: an increase in 
F max , a decrease in V max , and a characteristic increase 
in «/F max . 

An increase in [P ; ] reduces F max (Dawson et al., 1978; 
Nosek et al., 1987; Wilson et al., 1988; Potma et al., 
1995) without changing V max (Cooke etal., 1988), espe- 
cially at near in vivo temperature (Debold et al., 2004; 
Karatzaferi et al., 2008). The action of inorganic phos- 
phate on the cross-bridge cycle is more complicated 
than a two-state kinetic model can explain. It has been 
suggested that a partial power stroke can occur before 
the release of Pi from the cross-bridge (Dantzig et al., 
1992; Linari et al., 2010; Caremani et al., 2013). This 
could at least partially account for the depressive effect 
of Pi on F max . It has also been suggested that a myosin 
cross-bridge could detach from actin filament, even 
when the ATP hydrolysis products (P ; and ADP) are still 
attached to the bridge (Linari et al., 2010). This is a 
crucial feature used by Chin et al. (2006) to simulate 
and explain the unique effects of Pi that uncouples F max 
and V max . A simplified explanation is that P; reduces f&?v 
without affecting gxp P . This would imply a decrease 
in a/F max with increasing [P ; ], which is supported by 
Debold et al. (2004) but not Karatzaferi et al. (2008). 

In experiments where intracellular concentrations of 
ATP, ADP, and P; must be controlled, it is necessary to 
permeabilize the cell membrane. To preserve the struc- 
tural integrity of the "skinned" cells, the experiments are 
usually conducted at temperatures much lower than the 
body temperature. Readers must be aware that the results 
discussed above may be different at higher temperatures. 

Effects of muscle fatigue. Muscle fatigue is associated 
with reduced muscle performance, especially in power 
output, which in turn is determined by muscle force 
and velocity and the curvature of the force-velocity 
relationship. At the level of actomyosin interaction, 
muscle fatigue is likely a result of altered chemical 
environment caused by rapid use of ATP and a reduc- 
tion in intracellular [Ca 2+ ] (Fitts, 2008; Jones, 2010). 
The products of ATP hydrolysis include ADP, P i; and 
H + . Changes in intracellular [ATP], [ADP], [P ; ], [H + ], 
and [Ca' + ] associated with fatigue are likely direct modi- 
fiers of cross-bridge kinetics. To complicate the matter, 
the effects of some of these modifiers are temperature 
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dependent, for example, the effect of [H + ] (Knuth 
et al., 2006). To study mammalian muscle fatigue, ex- 
periments should be conducted at near body tempera- 
ture (Jones, 2010). The decrease in F max seen in fatigue 
is likely a consequence of high intracellular [P ; ] and 
[H + ], and low [Ca 2+ ]. P;, besides its direct effect on the 
cross-bridge cycle mentioned above, has the additional 
effect of reducing myofibrillar Ca 2+ sensitivity (Debold 
et al., 2006). High intracellular [H + ] is known to cause 
cross-bridges to accumulate in low force states (Seow 
and Ford, 1993; Debold et al., 2008), which could re- 
duce F max . Low intracellular pH is also a likely player in 
reducing V max during muscle fatigue. The myosin 
ATPase activity is lowered at low pH (Cooke et al., 
1988), perhaps because of lower rates of release of ADP 
from the myosin head (Debold et al., 2008), which has 
been found to reduce V max (Barany, 1967; Siemankowski 
et al., 1985; Elangovan et al., 2012) . It has been noticed 
that in the initial phase of muscle fatigue, small de- 
creases in F max could occur in the absence of apprecia- 
ble change in V max (Edman and Mattiazzi, 1981; Jones 
et al., 2006) . Because phosphocreatine is the energy re- 
serve (^20 mM) in muscle cells used first to replenish 
ATP in a process that produces inorganic phosphate, 
short-term fatigue is likely mainly caused by an increase 
in intracellular [P ; ] (Wilson et al., 1988). As mentioned 
above, high [P ; ] is able to reduce F max without affecting 
V max , making it a likely candidate for causing the initial 
decrease in F max seen in muscle fatigue. The cause for 
the reduction in a/F max (increase in force-velocity cur- 
vature) seen in muscle fatigue is less clear and more 
controversial (Fitts, 2008). However, the individual or 
combined effect of high [P ; ] and [H + ] is likely the cul- 
prit. Low [Ca 2+ ] appears to have no effect on a/F max 
(Langeron et al., 1999). 

Reduction in a/F max seen in muscle fatigue is a major 
factor contributing to the loss of muscle power (Jones 
et al., 2006). Because both elevated [P ; ] and [H + ] re- 
duce the cross-bridge transition rate from low to high 
force states (Cooke et al., 1988; Seow and Ford, 1993; 
Debold et al., 2004, 2008; Karatzaferi et al., 2008), the 
decrease in a/F max is likely a result of a reduction in fi^p 
(according to Eq. 12; and according to Fig. 3 B and Eq. 14, 
a reduction in leads to reduction in /y> P ). Jones 
(2010) has reached a similar conclusion that there must 
be a substantial decrease in the rate constant for attach- 
ment in the Huxley (1957) model to account for the 
observed decrease in power and increase in force-ve- 
locity curvature in fatigued muscle. 

Effects of myosin light chain phosphorylation in striated 
muscle. In submaximally activated striated muscle fibers, 
isometric twitch and tetanic force are substantially aug- 
mented by the phosphorylation of the regulatory myo- 
sin light chain (Persechini et al., 1985; Sweeney and 
Stall, 1986; Grange et al., 1995). One of the effects of 



light chain phosphorylation is that it pushes myosin 
heads toward the actin filaments, and this has been spec- 
ulated to facilitate the actomyosin interaction (Metzger 
et al., 1989; Levine et al., 1996; Yang et al., 1998). Analy- 
sis of data on isometric tension potentiation and redevel- 
opment (Metzger et al., 1989; Sweeney and Stall, 1990) 
suggests that the apparent rate of attachment (/app) is 
enhanced by light chain phosphorylation. Examination 
of the changes in force-velocity properties and power 
output caused by phosphorylation of the regulatory 
light chain reveals an increase in the a/F max ratio 
(Grange et al., 1995), consistent with the hypothesis 
that /app is enhanced by phosphorylation of the myosin 
light chain that increases the muscle power in part by 
reducing force-velocity curvature. 

Hill's hyperbola, no longer just descriptive 

Despite being regarded as purely empirical and lacking 
physiological insight soon after its inception three quar- 
ters of a century ago, the Hill hyperbola has continued 
to be the predominant descriptor of the mechanical 
performance of muscle. Although the connection be- 
tween the hyperbolic equation and cross-bridge cycling 
kinetics could have been made long ago, soon after AF. 
Huxley's publication of his 1957 model, a solid connec- 
tion was only possible after the state-of-the-art study on 
the mechanics of myosin motors by Piazzesi et al. (2007) 
at nanometer and piconewton resolution. For physiolo- 
gists and kinesiologists who use Hill's equation to ap- 
proximate force-velocity relationships, this connection 
provides molecular mechanistic insights to inform their 
interpretations of data. The connection also provides a 
new perspective for appreciating a seminal contribu- 
tion to muscle physiology that is of both historical and 
contemporary interest. 
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